Two-dimensional structure in a generic model of triangular proteins and protein 

trimers 
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Motivated by the diversity and complexity of two-dimensional crystals formed by triangular pro- 
teins and protein trimers, we have investigated the structures and phase behavior of hard-disk 
trimers. In order to mimic specific binding interactions, each trimer possesses on 'attractive' disk 
which can interact with similar disks on other trimers via an attractive square- well potential. At low 
density and low temperature, the fluid phase mainly consists of tetramers, pentamers, or hexamers. 
Hexamers provide the structural motif for a high-density, low-temperature periodic solid phase, but 
we also identify a metastable periodic structure based on a tetramer motif. At high density there is 
a transition between orientationally ordered and disordered solid phases. The connections between 
simulated structures and those of 2D protein crystals - as seen in electron microscopy - are briefly 
discussed. 

PACS numbers: 64.60.-i,64.70.Dv,05.10.Ln 



I. INTRODUCTION 

Two-dimensional (2D) materials present some fasci- 
nating challenges to condensed-matter theory, with even 
the most simple 2D systems harboring surprises. One of 
the most famous problems involves the precise descrip- 
tion of melting in 2D solids made up of hard, disk-like 
particles with short-range repulsive interactionsiS*^ 
Specifically, does the fluid undergo a weak first-order 
transition to the solid, or is there an intermediate hexatic 
phase linked by two continuous phase transitions? Re- 
lated avenues of research concern the existence of exotic 
phases in systems made up of more complex particles, 
such as (non)periodic solids of hard-disk dimers;^ pen- 
tamers and hexamers^ tetratic phases of hard squares? 
and hard rectangles^ and orientationally ordered solids 
of hard pentagons and heptagons^ The effects of addi- 
tional interactions on the phase behavior and dynamics 
of 2D systems are also of interest, as evidenced by recent 
studies on dipolar potentials in the context of magnetic 
colloids. 11 Such models provide an ideal testing ground 
for condensed-matter theories, and in some cases chal- 
lenge our most fundamental understanding of the prop- 
erties of matter. 

Despite their simplicity, 2D models can provide reliable 
descriptions of some real, and rather complex, experi- 
mental situations. For example, in a number of recent 
studies, 2D models have been employed to help inter- 
pret and understand the clustering and crystallization of 
proteins at interfaces. The conformations and interac- 
tions of proteins are central to biological activity, and 
ideally one would like to investigate these properties in 
vivo. Unfortunately, structural information is most com- 
monly obtained from X-ray diffraction studies on crys- 
tals. There is a class of proteins, however, that can be 
studied under conditions resembling those in vivo. Mem- 
brane proteins constitute a large class of molecules found 
within the lipid bilayers that constitute cell walls. They 
fulfill a variety of roles, such as controlling the selective 



transport of ions and molecules across cell membranes, 
or providing binding sites for other molecules on to the 
membrane. The structures of membrane proteins can be 
studied by deposition on to a surface, alongside lipids, 
to form either low surface-coverages or 2D crystals; the 
hydrophobic lipids help to mimic the interior of the mem- 
brane. Electron microscopy or atomic-force microscopy 
can then be used to image directly the clustering and 
packing of proteins at the solid-air interfaced In many 
cases, the ordering of proteins can be rationalized on 
the basis of their gross shapes (the way in which those 
shapes would 'tile the plane') and the presence of spe- 
cific binding interactions between domains on different 
molecules. For example, the surface structure of bacteri- 
orhodopsin (a transmembrane protein) is comprised of a 
close-packed array of trimers, each made up of monomers 
that resemble 120° sectors of a circle. Monte Carlo (MC) 
simulations of hard sectors - with an additional attrac- 
tive square-well potential to mimic specific binding in- 
teractions - yield insight on the self-assembly and subse- 
quent crystallization processes^ In another application, 
the ordering in 2D crystals of annexin V - another 'tri- 
angular' membrane protein - was reproduced in simula- 
tions of a hard-disk model decorated with an appropriate 
orientation-dependent potential to mimic the locations 
of the specific binding sites on the protein. Experimen- 
tally observed honeycomb and triangular structures were 
captured by the molecular model. These examples show 
that the basic physics of large-scale structural order in 2D 
protein crystals can be studied with simple models, and 
without resorting to atomistically detailed - and hence 
very expensive - computer simulations. 

There are a large number of proteins which are either 
inherently triangular in shape, or otherwise form trimeric 
structuresii2ii£ii§iiLi£ii2i2ii For example, rotavirus inner 
capsid protein V6 forms trimers resembling equililateral 
triangles, which pack in 2D crystals (space group p6)d 
Specific fragments of prion proteins found in vivo form 
trimeric units that crystallize in to a 2D structure (space 
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group p3) Finally, we highlight an example in which a 
membrane fusion protein (from the Semiliki forest virus) 
is seen to form pentagons of trimers, with the center of 
the pentagon raised slightly out of the plane. 18 Some sem- 
blance of local five-fold coordination can also been seen 
in TetA - a roughly triangular transporter protein - at 
moderate surface coverages^ 

Motivated by the diversity of 2D crystal structures ex- 
hibited by trimeric protein units, and also by the observa- 
tion of five-fold coordination^ we have investigated the 
structure and phase behavior of model trimeric molecules 
made up of hard disks. In order to mimic specific binding 
interactions, such as those that might give rise to local 
five-fold coordination, we focus on an equililatcral trian- 
gle of three hard disks at contact, in which one disk can 
interact with the corresponding disks on other molecules 
via a short-range attractive square- well potential. As we 
will show below, this raises the possibility of generating 
orientational order within simple close-packed structures, 
and also offers the opportunity of forming clusters at low 
surface coverages. Using MC simulations, we map out 
the phase diagram of the model system, and character- 
ize the structures of the low-density clustered fluid and 
high-density solids which are formed at low temperature. 
The remainder of the article is organized as follows. In 
Section^ we describe the molecular model, and summa- 
rize the simulation methods. The results are presented 
in Section ITTT1 and Section llVl concludes the paper. 

II. MODEL AND METHODS 

The molecular model consists of three hard disks, each 
of diameter a, fused at mutual contact to form an equi- 
lateral triangle. Two of the disks on each molecule are 
purely repulsive, and interact with all other disks in the 
system through the potential 

where r is the separation between the centers of two disks. 
The third disk on each molecule carries a central attrac- 
tive interaction site; these 'attractive' disks interact with 
each other via the potential 

!oo r < a 
-e cr < r < Act (2) 
r > Act 

where A > 1 controls the range of the attraction. This 
potential crudely mimics an effective attraction between 
vertices of the molecular triangles, which might arise 
through specific interactions (e.g. hydrogen bonding, 
disulfide bridges, effective solvophobic interactions). 

The parameter A will clearly have a crucial role to play 
in the thermodynamics of the system. If A >> 1 then 
one should anticipate a conventional phase diagram con- 
taining a vapor-liquid transition, and a fluid-solid tran- 
sition. The orientation of a trimer can be defined by 



a vector n joining the geometrical center of the trimer 
with the center of the attractive disk. It is unlikely that 
there would be any periodic orientational ordering of n 
in the solid phase; if two trimers can interact favorably 
irrespective of the mutual orientation, then on entropy 
grounds the orientations will be disordered. In the op- 
posite extreme, A > 1, the molecules will feel the ori- 
entation dependence of the net trimer-trimer potential, 
and ultimately we might expect the vapor-liquid transi- 
tion to disappear from the equilibrium phase diagram. 
Indeed, in a pure square-well hard-sphere fluid, conden- 
sation becomes metastable with respect to freezing when 
A < 1.25. 21 In the present case, an interaction range 
A < v3 guarantees that attractive sites must face each 
other directly in order to interact; when A > y3 it is 
possible for an attractive disk to be within interaction 
range of a trimer even if it approaches from 'behind'. 
With these comments in mind, we have chosen to study 
a system with A = 1.25. The ratio of Act to the (angle- 
averaged) diameter of the trimer is smaller than that in a 
pure square- well hard-sphere system with the same value 
of A, and assuming some sort of correspondence between 
two and three-dimensional systems, we do not anticipate 
there being a vapor-liquid transition in the equilibrium 
phase diagram. On the other hand, because the trimers 
have to attain quite specific mutual orientations in order 
to interact favorably (since A < y/S), we should expect 
to see some sort of non-trivial structure in fluid and solid 
phases at low temperatures. 

Systems of N = 120 trimers were studied using MC 
simulations either in the isothermal-isobaric (NpT) en- 
semble or the canonical (NAT) ensemble^ 2 . The simu- 
lation cell was rectangular with dimensions L x and L Vl 
and area A — L x L y . Each MC cycle consisted of one 
translational trial move and one orientational trial move 
for each of N randomly selected molecules. To help equi- 
librate dense phases, every fifth MC cycle included N 
trial moves in which a randomly selected trimer was ro- 
tated by ±120°. In NpT simulations of solid phases, L x 
and L y were varied independently; in NpT simulations 
of fluid phases, the simulation cell was constrained to 
be square. For most thermodynamic state points typi- 
cal equilibration runs consisted of ~ 10 5 MC cycles, but 
some points (close to phase transitions) required ~ 10 6 
MC cycles. Production runs were typically ~ 10 5 MC 
cycles. We define the following dimensionless units in 
terms of the square-well depth, e, and the hard-disk di- 
ameter, ct: number density p* = No /A; temperature 
T* = k B T/e; pressure p* = pa 2 /k B T. 



III. RESULTS 

The phase diagram of the model trimers in the density- 
temperature (p*-T*) plane is sketched in Fig. ^ Before 
detailing the determination of the phase boundaries, the 
characteristics of the different phases will be described. 
There are four distinct regions in the phase diagram. At 
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FIGURE 1 

FIG. 1: Phase diagram of the model trimer system in 
the density-temperature (p*-T*) plane: (solid points and 
solid lines) approximate fluid-solid phase boundaries, as- 
sumed to be first order; (open points and dashed lines) bound- 
aries between high-temperature unclustered states and low- 
temperature clustered states, as evidenced by maxima in the 
heat capacity along isobars; (dot-dashed line) close-packed 
density, p* p = 2/3V<$ ~ 0.3849. 



low density and high temperature, a normal fluid phase 
is in evidence (fluid I) . A typical simulation configuration 
is shown in Fig. |2Ja) ■ There is neither translational nor 
long-range orientational order in the system. 

At high density and high temperature, the stable solid 
phase (solid I) possesses an orientationally disordered 
structure (in the sense that n is disordered) with the 
trimers close-packed to form alternating rows displaced 
by a/2. Figure Hfb) shows both the lack of orientational 
order, and the registry between alternating rows. Notice 
the black bonds showing how the disks are connected 
within the trimers; we call this an l AB' 1 structure to de- 
note the alternating alignment of the rows. The close- 
packed rows resemble those formed by VP6*i^ although 
the registry between the rows is different. At the end 
of this section, we will briefly discuss the possibility of 
solids with other close-packed structures. 

At low temperature and low density we find a highly 
associated fluid (fluid II) , in which the attractive disks ag- 
gregate to form distinct clusters. A typical configuration 
is shown in Fig. Hfc), which exhibits a broad distribu- 
tion of cluster sizes. To identify clusters, we employ the 
obvious criterion that two trimers with attractive disks 
within interaction range belong to the same cluster. With 
this definition in mind, Fig. |2Ic) shows that, in general, 
the attractive disks within the clusters form close-packed 
motifs, rather than loose arrangements of disks on the cir- 
cumference of a ring. For clusters of three trimers there is 
no distinction, whereas for four or more trimers the close- 
packed arrangement is more favorable; in a ring, each disk 
would have two nearest neighbors, whereas close-packed 
motifs can accommodate more than two direct contacts. 



In Fig. |3| we show the probability distribution function 
of clusters containing n molecules, at different pressures 
along an isotherm with T* = 0.3. As the pressure and 
density increase, the distributions show peaks at progres- 
sively higher values of n. At the highest fluid-density 
shown - p* = 0.280, Fig. |3Je) - the most probable clus- 
ter size is n = 5. We had hoped that these clusters would 
adopt a pentagonal structure, but instead the attractive 
disks form 'Olympic rings' motifs, such as those shown 
in Fig. |2Jc). The maximum disk-disk separation in a 
perfect pentagon of disks is \/2(l — cos 108°)ct ~ 1.62a, 
which is longer than the range of the potential studied 
in this work. Hence, to minimize the energy, the cluster 
will contract to form a close-packed structure. Perhaps 
pentagonal clusters would be formed in a system with 
1.62 < A < a/3? (The upper limit means that there 
can be no other disks between two interacting attractive 
disks.) We did some test runs in the fluid phase with 
A = 1.7, but no pentagonal clusters were observed. If 
anything, fewer distinct clusters were in evidence as com- 
pared to A = 1.25, presumably because it is less crucial 
that the trimers attain a specific mutual orientation in 
order to interact. 

Upon compression of the low-temperature fluid we 
often encountered metastable structures, such as that 
shown in Fig. Hfd). This clearly shows a predominance 
of n = 6 clusters, with the attractive disks close packed 
to form a parallelogram motif, but the clusters are not 
yet fully packed in to a solid structure. This process is 
completed upon further compression, to form a p2 pe- 
riodic solid (solid II), a defective example of which is 
shown in Fig. [2[e). In simulations of the high-density 
solid II phase, the initial configuration consisted of the 
appropriate AB structure, but with n for each molecule 
chosen randomly from the three molecular arms; the ori- 
entational structure shown in Fig. Hje) develops spon- 
taneously. The cluster distribution for such a solid at 
temperature T* = 0.3 and density p* — 0.329 is shown 
in Fig- Etf) - The primary peak is at n = 6, but the pres- 
ence of defects - such as those shown in Fig. |2fe) - gives 
rise to smaller 'clusters' of attractive disks. 

The fluid-solid phase boundaries were located by moni- 
toring the equation of state p{p) along selected isotherms 
in NpT simulations. For each isotherm, two sets of sim- 
ulations were performed: a compression branch, start- 
ing from a low-density fluid configuration; and an ex- 
pansion branch, starting from the perfect solid structure 
corresponding to that found in the compression branch 
at high pressure. Portions of two representative ex- 
amples (T* = 0.3 and T* = 1) are shown in Fig. ^ 
Of course, the fluid equations of state extend to much 
lower densities, but these exhibit entirely conventional 
behavior and hence are not shown; in particular, there 
is no sign of a 'van der Waals' loop which would indi- 
cate a vapor-liquid phase transition. The main features 
of interest are the apparent discontinuities in the den- 
sity at what are assumed to be first-order phase tran- 
sitions (we will not open up the can of worms associ- 
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FIG. 2: (Color online) Configuration snapshots from NpT simulations: (a) normal fluid phase (fluid I) at T* = 2, p* — 2.5, 
p* = 0.259; (b) orientationally disordered A_B-solid phase (solid I) at T* — 2, p* = 12, p* = 0.345; (c) clustered fluid phase (fluid 
II) at T* = 0.25, p* = 0.75, p* = 0.222; (d) metastable state at T* = 0.25, p* = 2.6, p* = 0.290; (e) orientationally ordered 
AB-solid phase (solid II) at T* = 0.25, p* = 12, p* = 0.349; (f) metastable orientationally ordered AA-solid at T* = 0.25, 
p* = 20, p* = 0.356. In each case the attractive disks are colored dark gray (red online), the repulsive disks are colored light 
gray, and all disks are drawn with diameter la. 



ated with the precise nature of two-dimensional melting 
and freezingi* 2 * 3 *^) . In Fig. 0] we indicate distinct fluid 
and solid branches in the equations of state, a number 
of putative metastable states (as discussed above), and 
approximate tie-lines connecting the fluid and solid co- 
existence densities, obtained as follows. The fluid branch 
was fitted with a virial expansion containing terms up to 
p 5 , i.e., p/ksT — p + J2n=2 B n p n , while the solid branch 
was found to be fitted rather well by a simple van der 
Waals equation 23 of the form p/kT = ap/(l — bp) — cp 2 , 
which contains a free- volume term arising from repulsive 
interactions, and a mean-field term arising from the at- 
tractions. The coexistence densities were then estimated 
by extrapolating the fitted branches of the equation of 
state to a pressure half way between those in the highest- 
density stable fluid and the lowest-density stable solid; 



the metastable states were identified as those that did 
not fit on to either branch and/or for which the simula- 
tion configuration was clearly neither pure solid nor pure 
fluid, e.g. Fig. Hfd). Obviously this approach provides 
only very rough locations for the phase boundaries shown 
in Fig.^ but some general trends are nonetheless appar- 
ent. At very low temperatures, the coexistence densities 
decrease as the system is cooled, and the transition ap- 
pears to be getting weaker. At high temperatures (T > 1) 
the fluid coexistence density (ft* ~ 0.30) is very similar 
to the density at which the pure hard-disk fluid under- 
goes its transition, either to a hexatic or a solid (disk 
density p* = 0.899^ 'trimer' density p* = 0.300). The 
apparent trimer solid coexistence density (p* ~ 0.32) is 
significantly larger than the melting density of hard disks 
(disk density p* ~ 0.914^ 'trimer' density p* ~ 0.305). 
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FIG. 3: Cluster distributions for systems along the isotherm 
T* = 0.3: (a) p* = 0.5, p* = 0.188; (b) p* = 1, p* = 0.230; 
(c) p* = 1.5, p* = 0.252; (d) p* = 2, p* = 0.265; (e) p* = 2.5, 
p* = 0.280; (f) p* = 6, p* = 0.329. In (a)-(e) the system is 
fluid, whilst in (f) the system is solid (II). 
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FIG. 5: Configurational energy U (left) and excess heat ca- 
pacity Ca (right) as functions of reduced temperature T* at 
the close-packed density p* — 2/3%/3 ~ 0.3849: (circles) sim- 
ulation results; (lines) results derived from a Pade [5,5] fit (see 
text). 
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FIGURE 4 

FIG. 4: Equations of state along isotherms with T* = 0.3 
(solid symbols, solid lines), and T* = 1 (open symbols, dashed 
lines): (circles) state points from NpT simulations, with 
AB solid phases; (squares) state points from NpT simula- 
tions, with A A solid phases (T* = 1 only); (crosses) putative 
metastable state points; (triangles) approximate coexistence 
densities; (lines) fits to the fluid and solid branches (see text). 
The statistical errors in the NpT simulation points are smaller 
than the symbols. 

The final piece of the equilibrium phase diagram con- 
cerns the crossover from high-temperature orientation- 
ally disordered states to low-temperature states that pos- 
sess structural motifs arising from the clustering of the at- 
tractive disks. To delineate the boundary between these 
two regimes, we calculated the heat capacity appropri- 
ate to the statistical mechanical ensemble being sam- 
pled. In general we used NpT simulations to measure 
C p = (dH/dT) p - where H = U + pA is the enthalpy 



(minus the kinetic contribution) - as a function of tem- 
perature along an isobar. Since clustering must be ac- 
companied by a drop in the configurational energy, and 
enthalpy, a peak in C p would seem to be an obvious sig- 
nal of a crossover from unclustered to clustered states. In 
simulations we evaluated the usual fluctuation formula, 
C p = [(H 2 ) - {H) 2 ]/k B T 2 , and, as a check, differenti- 
ated an [n, n] Pade approximant fitted to the enthalpy as 
a function of T; 

= a + ai T + a 2 T 2 + ... + a n T n 

l + b 1 T + b 2 T 2 + ... + b n T n ' ( ' 

These two approaches yielded consistent results, and the 
peak in C p was easy to locate accurately. In general the 
peak height is less pronounced at high densities, mainly 
due to the fact that even in the high-temperature phase 
there must be some attractive disks within interaction 
range due to the confinement. Thus, the most difficult 
situation obtains at close packing of the trimers, p* p = 
2/3V3. In this case we studied a perfect close-packed AB 
solid, and carried out NAT MC simulations with ±120° 
rotations only. We show results for the configurational 
energy, U , and the excess constant-area heat capacity, 
C A = (dU/dT) A , in Fig.0 A [5,5] Pade fit provides a 
reliable description of the energy, and the corresponding 
results for Ca are consistent with those obtained via the 
fluctuation formula. 

In Fig. ^ we show the positions of the maxima in C p 
- and Ca at p* p = 2/3V3 - along with separate cubic 
fits to the points in the fluid and solid regions of the 
phase diagram. It appears that the two branches would 
meet up somewhere in the fluid-solid coexistence region. 
We stress that the boundaries indicated do not represent 
thermodynamic phase transitions; rather, they separate 
different regimes of trimer association. 



Finally, we briefly consider the possibility of the trimer 
system adopting other solid structures, such as the p2 AA 
structure shown in Fig. 0T), in which the close-packed 
(horizontal) rows are matched with the neighboring rows. 
In this case, the low-temperature, orientationally ordered 
solid exhibits rhombic cluster-motifs containing only four 
attractive disks. Out of those four disks, two are inter- 
acting with two other disks, and two are interacting with 
three other disks. Hence, the minimum configurational 
energy for an A A solid is — |e per trimer. In the AB 
structure, there are six attractive disks per parallelogram 
motif, of which two have two neighbors, two have three 
neighbors, and two have four neighbors, giving a min- 
imum energy of — le per trimer. Hence, on energetic 
grounds, we should expect the AB structure to be ther- 
modynamically favored. Even at high temperature, the 
AA structure appears to be less stable with respect to 
the AB structure. As an example, in Fig. we show 
an yL4-solid branch of the equation of state at T* = 1, 
alongside the AB-solid branch. For a given pressure, the 
AB solid has the higher density which makes this state 
at least mechanically stable with respect to AA. Indeed, 
we only ever observed the fluid spontaneously freezing in 
to an AB structure. Although we have not performed 
free-energy calculations, it would be very surprising if an 
cntropic effect could compensate for the relative energetic 
and mechanical stability of the AB phase with respect to 
the A A phase. 

Another possible close-packed structure is illustrated 
in Fig. inia), without any indication of the attractive 
disks. This structure resembles that adopted by 2D crys- 
tals of TetAjiS although we never saw this packing struc- 
ture emerge from our simulations. As far as our model is 
concerned, the absence of this structure at low temper- 
ature is easy to understand. In Figs. Efb) and HJc) we 
illustrate mirror images of the most obvious periodic ar- 
rangement of the attractive disks (space group p3). The 
energy per trimer is only — le, and so this is not compet- 
itive with the AB structure that is seen to emerge spon- 
taneously in our simulations. Free-energy calculations 
would be of interest, particularly at high temperatures 
where entropy is everything! 



IV. DISCUSSION 

In this article we have described the structure and 
phase behavior of a generic model of trimeric molecules, 
largely motivated by recent experimental 2D microscopy 
studies of clustering and crystallization in triangular pro- 
teins and protein trimers. The molecular model consists 
of a triangle of hard disks, with one of the disks partic- 
ipating in attractive square-well interactions with simi- 
lar disks on other trimers. The range of the square- well 
potential, Act, was 1.25 times the disk diameter. This 
system crudely mimics the general shape and specific in- 
teractions of a wide range of proteins. The model system 
exhibits fluid and solid phases which, at low tempera- 
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FIG. 6: (Color online) Illustrations of an alternative close- 
packed structure: (a) without an assignment of attractive 
disks; (b) and (c) mirror images of a possible structural motif 
for a periodic arrangement of attractive disks. The attractive 
disks are colored dark gray (red online), the repulsive disks 
are colored light gray, and all disks are drawn with diameter 
la. 



tures, possess interesting structural motifs arising from 
the clustering of the 'attractive' disks. 

In the fluid, a distribution of clusters is in evidence, in- 
cluding tetramers, pentamers, and hexamers (of trimers). 
In the pentamers and hexamers, the attractive disks 
close-pack to form 'Olympic rings' and parallelogram 
shapes, respectively. We had hoped to find more open 
pentagonal clusters of trimers, such as those reported in 
Ref. [l8l To investigate the formation of such clusters 
further, it might be interesting to study a system of hard 
isosceles triangles with the unique angle equal to 72°, 
and a short-range attraction operating between the cor- 
responding vertices. 

In the low-temperature solid, the basic structural motif 
consists of clusters of six molecules, with the attractive 
disks close-packed to form a parallelogram. A metastable 
solid possessing a motif made up of four molecules was 
also identified. The fundamental difference between the 
two situations is the registry between neighboring close- 
packed rows of trimers (AB versus AA). Even at high 
temperatures, the orientationally disordered AB solid is 
at least mechanically stable with respect to the AA solid. 
We identified a third structure based on hexagonal close 
packing, but this structure is not competitive either, at 
least in terms of energy. It would be worth performing 
free-energy calculations to study these issues further. 

Finally, it is worth commenting that a diverse range of 
2D structures can be generated from very simple molec- 
ular models. Fully atomistic calculations of 2D protein 
structures are expensive, and, it could be argued, yield 
little insight on the fundamental physics behind cluster- 
ing and crystallization. As has been shown in a variety 
of cases, including the present study, the process of de- 
veloping and studying simple models of complex systems 
can yield some surprising results. 
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